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In models with large extra dimensions all gauge singlet fields can in principle propagate in the extra 
dimensional space. We have investigated possible constraints on majoron models of neutrino masses 
in which the majorons propagate in extra dimensions. It is found that astrophysical constraints 
from supernovae are many orders of magnitude stronger than previous accelerator bounds. Our 
findings suggest that unnatural types of the "see-saw" mechanism for neutrino masses are unlikely 
to occur in nature, even in the presence of extra dimensions. 
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I. INTRODUCTION 

Over the past few years there has been an enormous in- 
terest in models with extra dimensions beyond the 3+1 
dimensions of the standard model. Such extra dimen- 
sions can be large if the gauge fields of the standard 
model are constrained to a 3+1 dimensional brane in 
the higher dimensional bulk space. This idea is particu- 
larly interesting because it may be connected with a low 
quantum gravity scale, perhaps only slightly above the 
scale of electroweak symmetry breaking. In Refs. |l|. |[ ||] 
such a model was proposed in which the standard model 
brane is embedded in a bulk space of dimension S which 
is toroidally compactificd. Because gravity is allowed to 
propagate in the full bulk space it will look weaker to an 
observer confined to the brane. By Gauss' law the effec- 
tive four-dimensional Planck scale, Mp, can be related 
to the true A+5 dimensional energy scale of gravity, M , 
by 



Mp = {2nR) s M 



8+2 



(1) 



where R is the common radius of all the extra dimen- 
sions. If R is sufficiently large M can be as small as the 
electroweak scale. Effectively gravity is weak because the 
field lines leak into the extra dimensions. 

One of the most interesting features of these models is 
the presence of a tower of new modes for particles propa- 
gating in the bulk. The momentum of particles propagat- 
ing in extra dimensions is discretised because the extra 
dimensions are compact. To an observer on the brane 
this effectively looks like a tower of new states (Kaluza- 
Klein) for each bulk particle, with a 4D mass related to 
the extra dimensional momentum. The energy spacing 
between these states is then in general of order R^ 1 which 
can be very small. 

For the graviton this property has been used to put 
tight constraints on the possible radii of the extra di- 
mensions. However, any particle which is a singlet in the 



standard model could in principle propagate in the bulk. 
Examples of this are sterile neutrinos and axions. 

Some authors, for example, have investigated the pos- 
sibility that right handed neutrinos might propagate in 
the extra dimensions and that their resultant suppressed 
couplings to the usual left handed neutrinos could ac- 
count for the low scale of neutrino masses H, H , H , [?], ||, 
§, 0, fO). Other authors @ || 0, || p|, ]LQ |§ have 
studied the constraints on such models due to experiment 
and observation. An alternate approach |is| , correspond- 
ing to the conventional see-saw mechanism, has also been 
investigated by some authors. In this approach a Higgs 
singlet, which carries no standard model gauge quantum 
numbers, is allowed to propagate in the extra dimensions 
and give mass to right handed neutrinos which, for sim- 
plicity, are assumed to live on the brane. 

Recently in Ref. |Q a model was studied in which 
the lepton number is spontaneously broken and the asso- 
ciated Goldstone boson (the "Majoron") is present (see 
also j2l) ) . This model was used to compute the decay rate 
for the intermediate vector boson Z to two neutrinos and 
the Majoron (denoted by J) or one of its "Kaluza-Klein" 
excitatio ns |20l . Using the accurately known value of the 
Z width pa bounds on the dimensionless coupling of the 
two neutrinos to the Majoron were obtained allowing for 
any number of extra dimensions and any intrinsic mass 
scale (see also 23 ) . The related neutrinoless double beta 
decay process n + n — ► p + p + er + e~ + (J) has already 
been treated in a model of the present type 19 . A de- 
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tailed discussion of supernova constraints in the 3 + 1 
dimensional theory has very recently been given in (24). 

In this paper we use this model ]2C| to compute pro- 
cesses of astrophysical interest. We first see that when 
compact extra dimensions are present new processes be- 
come relevant and can heavily affect supernovae dynam- 
ics. We then show that supernovae constraints on the 
dimensionless couplings are by many order of magnitude 
more stringent than the accelerators bounds. Our find- 
ings seem to discourage unnatural see-saw models of neu- 
trino masses still allowed by accelerator bounds. 

We summarize in section II the Majoron model ex- 
tended to the extra-dimensional brane-world . In section 
III we briefly review the accelerator constraints. The Su- 
pernovae constraints are presented in section IV. In the 
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concluding section V we also discuss the possibility of the 
extra-dimensional Majorons to be source of dark matter. 



II. A MAJORON MODEL IN EXTRA 
DIMENSIONS 

In the original Majoron model of Chikashige, Mo- 
hapatra and Peccei |25) the lepton number associated 
with massive Majorana neutrinos is spontaneously bro- 
ken. Here, the notations of Ref. |26[ and Ref. |^(| for 
the 3+1 and the extra dimensional case will be followed, 
respectively. In addition to the usual Higgs doublet 



1 = 



(2) 



which has lepton number I equal to zero, the model con- 
tains an electrically neutral complex singlet field 



I 



(3) 



It is required that the Higgs potential constructed from <!> 
and </> conserves lepton number. The vacuum expectation 
values are: 



/J.0\ 



<**> , 

/J.0t 



(4) 



where Gf is the Fermi constant and ($) (whose non- 
zero value violates lepton number) sets a new scale in 
the theory. 

The three light physical two component neutrinos 
z^i , 1^2 , vz acquire Majorana masses mi, mi, m-$ which are 
of order e 2 M.H with 



O 



( v 



(5) 



According to the a standard "see-saw mechanism" ]27j] 
T>/\ represents the "Dirac-type" coupling constants for 
the bare light neutrinos while the 3x3 matrix A4h/{&) 
represents the Majorana type coupling constants for the 
bare heavy (or "right handed" ) neutrinos. Assuming the 
heavy scale Mh to be substantially larger than the en- 
ergy in play in the following we focus our attention only 
on the light neutrinos. The coupling of the Majoron J, 
identified as J = Im$, to the physical neutrino fields 
v \i v ii v z in 3 + 1 dimensions is: 



v a o 2 g ab v h 



(6) 



a,b=l 



It turns out |26| that the coupling constants have the 
expansion: 



gab 



-^rn a S ab 



(7) 



where the leading term is diagonal in generation space. 
It is convenient to express this leading term using four 
component ordinary Dirac spinors 



(8) 



in a 75 diagonal representation of the Dirac matrices the 
relevant Lagrangian term reads: 



J 



2($) 



i 



^ 5 ; , 7 1 

: Wa +4>a — 



75 



.(9) 

Here C is the charge conjugation matrix of the Dirac 
theory. 

The generalization of the present model to the case 
where the field $ propagates in 8 extra dimensions has 
been carried out in detail in Ref. [pofl . These extra di- 
mensions, denoted as yi with i = 1, . . .8, are assumed to 
be toroidally compactified with a radius Ri . For simplic- 
ity we assume that all the radii Ri are equal to the same 
value R. <& {x, y) continues to carry the "engineering di- 
mension" one as it would in 3 + 1 dimensional space-time 
and via a Fourier expansion with respect to the compact- 
ified coordinates we have: 

$(z,y)=Norm ^ * ni ,...,„, (x) e 4(™+-> , (10) 



where Norm = \2ttM R] 5 and M represents the intrinsic 
scale of the new theory. 

Each general Kaluza-Klein field receives a mass 
squared increment 



Ami 



1 

R? 



(11) 



The zero-mass Majoron Jo,o....,o( :c ) corresponds to the 
previously studied 3 + 1 massless Majoron. The fields 



.(x) are expected to receive a substantial in- 

and 



Re $ ni 

crement from the pure Higgs sector of the theory 
will be neglected in the following. 

The intrinsic scale M and the compactification ra- 
dius R can be related to each other when assuming in 
the "brane" model the graviton to propagate in the full 
(3 + 8) + 1 dimensional space-time. Then the ordinary 
form of Newtons' gravitation law is only an approxima- 
tion valid at distances much greater than R. The Newto- 
nian gravitational constant (inverse square of the Planck 
mass M P ) is obtained §], |, | as a phenomenological 
parameter from 



M P 
~M 



(2-kMRJ 



1 



(Norm) 



(12) 



Considering Mp as an experimental input (and approx- 
imating R\ = i?2 = • • • = Rs), shows via ( |12| ) that M is 
the only free parameter introduced to describe the extra 
dimensional aspect of the present simple theory when 8 
is fixed. 
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We expect Mh/{$) to be very roughly of the order 
unity and (<&) of the order of M. Finally the Yukawa 
interactions of the Majoron and its Kaluza-Klein excita- 
tions with the light neutrinos are described by (c.f. (ph): 

3 

£j — ^ ^ ^ ^ ^aa;ni,...,nj<^ni,...,n; ^ 
a=l ni,...,ns 

x^C^i^^ + ^^CC) ,(13) 
to leading order in the neutrino masses, m a with 



in,, 



M 



Qaa,' 



9a 



2($) M P ' 



(14) 



The vacuum expectation value (<&) (unless unnatural 
fine tuning is considered) is very roughly of the or- 
der of M, leading g aa to be naturally of the order of 
m a /Mp « 10~ 28 regardless of the number of extra di- 
mensions Clearly the exact value of this universal 
coupling crucially depends on the unknown dynamics be- 
hind the spontaneous breaking of the lepton number p0[ . 



III. REVIEW OF CONSTRAINTS FROM 
ACCELERATORS 

In |2(| it has been shown that the accurately known 
value of the Z width can provide information about the 
coupling of two neutrinos to the Majoron. Both the 3 + 1 
dimensional case and the case in which one adopts a 
"brane" world picture with the Majoron free to expe- 
rience the extra dimensions have been studied. Bounds 
on the dimensionless coupling constants were obtained, 
allowing for any number of extra dimensions and any in- 
trinsic mass scale. If the uncertainty of the Z's invisible 
width 1.7 x 1CP 3 GeV is roughly taken as an indication 
of the maximum allowed value for the total width into 
a Majoron and two neutrinos the following bounds on 
| g aa | were obtained in [§0| for M s = 10 4 GeV: \ g aa \< 
3.4 x 1CT 12 , | g aa |< 2.3 x 1CT 10 , | g aa |< 1.5 x 1(T 8 for 
5 = 2,3, 4, respectively. These are much stronger bounds 
than the one obtained for the model in 3 + 1 space-time 
dimensions which is | g aa |< 0.11. 

If a technically natural see-saw model is adopted, the 
predicted coupling constants are far below these upper 
bounds. In addition, for this natural model, the effect 
of extra dimensions is to decrease the predicted partial 
Z width, the increase due to many Kaluza-Klein excita- 
tions being compensated by the decrease of their common 
coupling constant. 

We shall see in the following that constraints from su- 
pernovae are much more stringent than the ones above. 



IV. CONSTRAINTS FROM SUPERNOVAE 

Supernova cooling — The proto-neutron stars created 
by core-collapse supernovae are born with core tem- 



peratures of 30-50 MeV. The main cooling mechanism 
for these stars is thermal surface emission of neutri- 
nos on a timescale of a few seconds. The total en- 
ergy emitted is the order a few 10 53 erg, with roughly 
equal amounts in all flavours of neutrinos and antineu- 
trinos. This emission has been observed from SN1987A 
by Kamiokande @], 1MB @ and Baksan all ob- 
serving v e events. The results from SN1987A are com- 
patible with theoretical supernova models and therefore 
put a constraint on any non-standard cooling mecha- 
nism carrying away too much energy, and since ma- 
jorons from the KK-tower will be produced in the su- 
pernova and carry away energy 1987A data can be used 
to constrain g aa . This has been done several times in 
the literature for the 3+1 dimensional majoron mod- 

| §, g [§. For 



els 



32 



3] £4|, |3J, |3§ j37j £8|, 

these models, only a relatively small band in parame- 
ter space is excluded for the following reason: For large 
g aa the majorons are tightly coupled inside the neutron 
star, and only escape via surface emission, like neutri- 
nos. Therefore they cannot carry away most of the en- 
ergy and supernovae do not yield significant constraints. 
On the other hand, for small g aa the majorons do es- 
cape freely, but they are not produced in significant 
numbers. The end result is that a band of roughly 
3 x 10~ 7 < g aa < 2 x 10~ 5 is excluded ||l). 

For the 3 + 1 + 5 models the situation is different be- 
cause each KK-mode is very weakly coupled. Therefore 
we are always in the limit where majorons escape freely 
once they are produced, and do not need to worry about 
possible surface emission. 

A fairly robust constraint on such "bulk emission" is 
the one proposed by Raffclt Q, that the emissivity of 
the neutron star medium should be 



e < 10 ergg s 



(15) 



In the 3+1 dimensional case where majorons are 
strictly massless the neutrino pair annihilation vv — > J 
is not allowed kinematically and the most important pro- 
cesses are vv — > J J and v — > D J. However, in the present 
case, most of the emitted majorons have mass of order T, 
and the simple pair-annihilation is by far the most im- 
portant process. The squared and spin-averaged matrix 
element for this process (per single Kaluza-Klein excita- 
tion) is 



1 



(16) 



The emissivity per volume of the medium is then 



Q(mj) = j dpidp 2 dpj(2Tr) 4 f v {pi)fv{p2) 

x\A\ 2 S 4 ( Pl +p 2 -pj){E 1 +E 2 ), (17) 

where dp = d 3 p/(2ir) 3 2E and /„ are the thermal distri- 
butions of the incoming neutrinos, f u = e~ p "/ T . Doing 
the integral gives 

Q(mj) = I ^-3 ff iT 5 x 4 K 2 (x) , x ^ mj/T , (18) 
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where K.2{x) is a Bessel function. However, we also need 
to sum over the KK-tower in order to obtain the total 
emissivity 



T 6/2 



M 2 T 5+«5 ^ 



Q = W2)W M 2+s & J dXXK ^ < 19 > 

= ^'^^^rs^ + m (20) 

Using the rough equality Q = e(p), (p) ~ 4x 10 14 gem -3 , 
this can be translated into a bound on g aa and M which 
is 



9aa 



^AM TeV i 30 



S)/2 



with the following values for X 



X 



1.0 x 10~ 21 

2.1 x 10~ 17 
4.6 x 10 
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for S = 2 
for 5 = 4 
for 5 = 6 



(21) 



(22) 



In the above, T30 = T/(30 MeV). In all figures we use 
T 30 = 1 . The bounds are shown in Fig. 1 as a function of 
M. However, from considerations of graviton emission, 
strong bounds on M already exist (44 , 45], p6| |47|, p), p9|. 
The most recent, and most restrictive, limits are those 
from Ref. j49|. These are shown as the thicker lines in 
the left part of the figure. Values of M in this range are 
already excluded from graviton emission arguments [ fi"9"| . 

Old neutron star excess heating — When consider- 
ing graviton emission a much stronger bound on M can 
be obtained by considering the decays of KK-gravitons 
produced in supernovae. Gravitons have a significant 
branching ratio into photons and the decays there- 
fore produce gamma rays. All cosmological supernovae 
have therefore by the present produced a diffuse cosmic 
gamma background. Comparing this with observations 
of the diffuse gamma background by the EGRET instru- 
ment |50|| yields a bound on M significantly stronger than 
found from the supernova cooling argument Jt8[ . More 
interestingly, the tightest constraint comes from observa- 
tions of old, isolated neutron stars. 

Most of the KK-modes emitted from the proto-neutron 
star have masses of order 3T and therefore also relatively 
low velocities. This again means that a large fraction 
(roughly one half) of the KK-modes have velocities lower 
than the escape velocity, and that neutron stars retain 
a halo of KK-modes with a typical radius of 2-3 Rns- 
When these gravitons decay they heat the neutron star 
and can lead to excess surface emission from old neutron 
stars. This argument applies equally to KK-gravitons 
and majorons. For gravitons it was used in Ref. jl9| to 
put an extremely strong limit on M, M > 1600 TeV 
(S = 2), 70 TeV (6 = 3). 

The KK-majorons only decay into neutrinos and not 
photons. However, the typical energy of the emitted neu- 
trinos is 50-100 MeV, and the neutron star is not trans- 
parent to neutrinos of such high energy. Therefore the 
neutrinos hitting the neutron star surface will heat it just 



like photons do. The surface luminosity of the neutron 
star at late times should therefore reach a constant level 
corresponding to the energy deposited by decay neutri- 
nos. This luminosity is given by 



LnS = fjFjETOT^j) 



R 



NS 



R 



(23) 



halo 



where fj is the fraction of the total supernova energy 
emitted in majorons, Fj is the fraction of the emitted 
majorons remaining bound to the neutron star, taken to 
be 1/2, Etot is the total SN energy, taken to be 3 x 10 53 
erg, (r) is the average majoron decay rate, and Rhaio 
is the typical radius of the majoron halo, taken to be 
2 Rns- The decay rate for non-relativistic majorons is 
given by 



1 

64^ 



(24) 



fj is a function of M and g aai and can be found from 
the above cooling bound. The cooling bound, to a good 
approximation, corresponds to half the total energy be- 
ing emitted in majorons (fj ~ 1/2). For gravitons 
the strongest bound comes from the neutron star PSR 
J0953+0755 |5l] ||, ||, which is the oldest neutron star 
for which the thermal surface temperature has been mea- 
sured. Its total surface luminosity is estimated to be 



10~ 5 L. 



For majorons this neutron star also 



(25) 



(26) 



yields a strong constraint on g aa which is roughly 

gaa * IVI TcV 1 30 ' 

with the following values for Y 

( 3.3 x 10~ 22 for 8 = 2 
Y = I 4.8 x 10~ 20 for S = 4 
I 7.0 x 10" 18 for 6 = 6 

In all cases this bound is significantly stronger than the 
cooling bound, just as it is for gravitons. The bounds 
are summarized in Fig. ^. However, there is a limit to 
the applicability of the bound. The age of the neutron 
star PSR J0953+0755 is estimated to be r = 1.7 x 10 7 
yr |5^]. If the majorons decay faster than this, they will 
have vanished already and cannot act as a heating source. 
From Eq. (J24F) one finds a lifetime of 



T=1.3xHr 21 <?- 2 



100 MeV \ 



TO J 



(27) 



giving a rough upper limit on g aa of 1.5 x 10~ 18 , above 
which the limit of Eq. ( p5| ) does not apply. 

Another possible problem is that majorons could be 
reabsorbed when they pass through the neutron star on 
a timescale much shorter than the age of the neutron 
star JB3]. The most relevant process for reabsorption is 
inverse bremsstrahlung, JNN — > TV TV which is induced 
via the electroweak interactions |2j| . Since the majoron is 
a pseudo-scalar one can estimate the rate for this process 
in the same way as for axions |43| ]. The result is that 
reabsorption only happens on a timescale much longer 
than the lifetime of the neutron star PSR J0953+0755. 
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FIG. 1: Upper bounds on g aa from the supernova cooling 
bound, Eq. (£lj), for various 8 and M. The bottom curve is 
for 8 = 2, the middle for 8 = 4, and the top for 8 = 6. The 
thick lines at low M correspond to the excluded region of M 
from graviton effects. 
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FIG. 2: Upper bounds on g aa from the neutron star heating 
bound, Eq. (^), for various 8 and M. The bottom curve is 
for 8 = 2, the middle for 8 = 4, and the top for 8 = 6. The 
thick lines at low M correspond to the excluded region of M 
from graviton effects. The horizontal line corresponds to the 
upper limit of applicability of the neutron star heating limit. 



V. DISCUSSION 
A. Majorons as dark matter? 

By the same neutrino pair annihilation process as in 
supernovae, majorons are also produced in the early uni- 
verse. This means that in principle one might obtain 
non-trivial bounds on g aa from considering cosmological 
production of majorons. For gravitons such considera- 
tions lead to extremely stringent bounds on AI and the 
maximum temperature, Trh, of the radiation dominated 
epoch after inflation [^5| ^TJ, |57], If the fundamental 



scale M is close to 1 TeV then the reheating temper- 
ature needs to be extremely low, typically in the MeV 
regime. However, there is a lower limit to how low Trh 
can be without disturbing big bang nucleosynthesis. De- 
tailed calculations have shown that this limit is roughly 
Trh > 0.7 MeV (H, & Thc number density of ma- 
jorons in the universe can be found by solving the inte- 
grated Boltzmann equation 

hj = -3Hnj + J^m^TK^mj/T), (28) 

which applies to a single majoron mode with mass mj. 
nj is the number density of the single majoron mode J 
and H is the Hubble parameter. By summing over all 
KK-modcs of the majoron in thc same way as it was 
done for gravitons in Refs. |H| [56|, one finds a present 
day density of 

xg 2 aa / dzz 5 - 1 / dqq 3 K x {q) (29) 

JO J z 

Requiring that this density is smaller than the critical 
density, p c = 8.1 x I0~ i7 h 2 GeV 4 , then yields an upper 
bound on g aa as a function of M and Trh- 

However, there is again a limit to the applicability of 
the bound because for large g aa the majorons will have 
decayed by the present. For the decay rate given in 
Eq. ( p4| ) and the requirement that r > 10 10 yr one finds 

9aa < 3 X 10- 19 T^ MoV , (30) 

assuming that the typical majoron mass is ~ ?>Trh, a 
fairly good approximation. In Fig. 3 we show the con- 
tours of g aa corresponding to critical density. Also shown 
is the lower bound on M as a function of Trh from con- 
sidering the decay of gravitons produced in the early 
universe. From this argument anything to the left of 
the thick grey lines in the figure is excluded. The lim- 
its used are the ones from Ref. |36) which are the most 
restrictive cosmological limits, combined with the limits 
from Ref. which for low Trh can be more restric- 
tive. Finally, we also show the upper bound on g aa from 
the above equation. In the case where the line from the 
graviton bound is to the right of the decay lifetime bound, 
majorons cannot possibly contribute to critical density. 
This is seen to be the case for both 8 = 2 and 6 = 6, and 
indeed for all values of S from 2 to 6. Therefore the con- 
clusion is that no non-trivial bound on g aa comes from 
cosmology, and further that majorons cannot contribute 
the dark matter of thc universe [[H] . Already, gravitons 
are excluded as a dark matter candidate because if they 
were to contribute critical density the photons produced 
by their decay would have been clearly visible. 
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FIG. 3: Contours of \o%{g aa ) which corresponds to p.j — p c . 
The upper panel is for 5 = 2 and the lower for 5 = 6. The 
thick line to the right is the lower bound on M coming from 
considering graviton emission. The thick vertical line to the 
left correponds to the maximum g aa for which majorons live 
longer than the age of the universe. 



B. Conclusions 



In the present paper we have shown that supernova 
constraints on brane-world scenarios for neutrino masses 
are many orders of magnitude stronger than the acceler- 
ator bound [^0| . Even so the constraints found are some- 
what weaker than what is naturally expected in see-saw 
models of neutrino mass. 

We have also shown that bulk majorons cannot act as 
the dark matter in the universe, at least not within the 
present scenario with equal radii of all the extra dimen- 
sions. 

Finally, our findings suggest that unnatural types of 
the "see-saw" mechanism for neutrino masses are un- 
likely to occur in nature, even in the presence of extra 
dimensions. 
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